Threshold feedback control for a collective flashing ratchet: threshold dependence 
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We study the threshold control protocol for a collective flashing ratchet. In particular, we analyze 
the dependence of the current on the values of the thresholds. We have found analytical expressions 
for the small threshold dependence both for the few and for the many particle case. For few particles 
the current is a decreasing function of the thresholds, thus, the maximum current is reached for zero 
thresholds. In contrast, for many particles the optimal thresholds have a nonzero finite value. We 
have numerically checked the relation that allows to obtain the optimal thresholds for an infinite 
number of particles from the optimal period of the periodic protocol. These optimal thresholds for 
an infinite number of particles give good results for many particles. In addition, they also give good 
results for few particles due to the smooth dependence of the current up to these threshold values. 

PACS numbers: 05.40.-a, 02.30.Yy 



I. INTRODUCTION 



Ratchets or Brownian motors are rectifiers of thermal fluctuations. This rectification is usually achieved through 
the introduction of an external deterministic or stochastic perturbation in a system that is or becomes asymmetric 
under spatial inversion lj. Over the last years ratchets have been studied due to their theoretical and experimental 
relevance. From a practical point of view the ratchet effect has many potential applications in biology, condensed 
matter and nanotechnology QQ- 

Ratchets can be viewed as controllers that act on stochastic systems with the aim of inducing directed motion 
through the rectification of the fluctuations. In particular, flashing ratchets are thermal fluctuation rectifiers based on 
switching on and off a periodic potential Several studies deal with the problem of optimizing the particle current 

[f| or the efficiency [f| in these systems. However, they all consider only open-loop controllers (as that obtained with a 
periodic or random switching) . Recently, feedback control protocols have been introduced in the context of collective 
ratchets 0, • In the feedback control protocols the action of the controller depends on the state of the system. This 
feedback control, or closed-loop control, can be implemented in systems where the state of the system is monitored 
(as occurs in some experimental setups with colloidal particles 0). 

In this paper we study one of these closed-loop controls, the threshold control, previously introduced in Ref. 
The structure of the paper is as follows. In the next Section we present the mathematical model of the collective 
flashing ratchet with the threshold control protocol and we discuss briefly other protocols that have been studied in 
recent articles. Later, in Sec. IIIII we analyze the dependence of the average center-of-mass velocity on the thresholds, 
obtaining analytical approximated expressions that are compared with the numerical results. In Subsec. IIII Al we 
study the small thresholds dependence (distinguishing the many particles case and the few particles case), while in 
Subsec. Till Bl we discuss the dependence of the average center-of-mass velocity for any thresholds and any number of 
particles. Finally, we present our conclusions in Sec. IIVI 
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II. THE MODEL 

We consider N Brownian particles with positions Xi(t) at temperature T within a ratchet potential V{x), and whose 
dynamics is described by the overdamped Langevin equations 

j±i(t) = a(t)F( Xi (t)) + Zi(t); i = 1, . . . , N, (1) 

with 7 the friction coefficient (related to the diffusion coefficient D through Einstein's relation D = ksT/^) and 
Gaussian white noises of zero mean satisfying the fluctuation-dissipation relation (£,i(t)£,j(t')} — 2^kBT5ij5(t — t 1 ). 
The force is given by F(x) ~ — V'(x) and a is a control parameter that can take only two possible values, a = 
(potential 'off') or a — 1 (potential 'on'). 

Several control strategies have been studied in order to maximize the particle current in this system. The optimal 
periodic switching 0, consists on switching the potential on during a time period T on and switching it off during 
T g. Note that it is an open-loop control protocol and therefore the results are independent of N. This protocol is the 
periodic flashing ratchet, that has been widely studied both theoretically and experimentally 0,0- The maximization 
of the center- of-mass instant velocity protocol has been introduced and studied in Ref. |3- It consist on switching 
the potential on only if the net force would be positive. Therefore, it is a closed-loop control protocol, because it 
needs information about the state of the system in order to operate. This is the best strategy for a single particle. 
However, for a large number of particles the system gets trapped with the potential 'on' or 'off' and then the average 
steady state current tends to zero as N increases ■ Another closed- loop control protocol, the threshold control, was 
proposed in 8] to avoid this effect. In this paper we analyze it further. 

The net force per particle is 



1 N 

/w = ^E F ^))- ( 2 ) 



On the other hand, given the state of the system Xi(t), a good estimator for the time derivative of f(t) can be obtained 
using Langevin equation £Q) and Ito calculus (see Ref. [3), 

i i 

The maximization of center-of-mass instant velocity protocol has a(t) = 0(/(t)), with O the Heaviside function 
[9(x) = 1 if x > 0, else 6(x) = 0]. In contrast, the threshold control policy has two thresholds u n > and u g < 
which induce earlier switchings that permit to avoid the trapping. When f(t) decreases below w on we switch off the 
potential, although the net force is still positive, in order to avoid the trapping. Analogously the potential is switched 
on if the net force per particle increases above u ff, so we induce the flipping of the system before f(t) is positive. 
Therefore, the threshold control is given by 



f 



if f(t) > 



_ i 1 if u off < f(t) < Uon and f cxp (t) > 0, ,^ 
| if Uoff < f(t) < u on and / exp (*) < 0, 
[0 if/(t)<u off . 

This scheme removes the long decaying tails in the evolution of the net force preventing the trapping. Note that this 
protocol and the maximization of the center-of-mass instant velocity protocol are feedback controls or closed-loop 
controls. The threshold control protocol in the zero thresholds limit gives the maximization of the center-of-mass 
instant velocity protocol. 

III. THRESHOLD CONTROL STRATEGY 
A. Small thresholds 



In this subsection we analyze the threshold control strategy improving and extending the analytic expressions found 
for the maximization of the center-of-mass instant velocity protocol • 
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1. Many particles: quasideterministic approximation 

For many particles (large N) the net force has a quasideterministic behavior. It can be described in terms of two 
contributions, a deterministic contribution /°° (given by the behavior for an infinite number of particles) plus a small 
stochastic contribution 

/(*) = f°°(t) + fluctuations. (5) 

This approximate description has proven to be fruitful in order to understand the behavior of these ratchets in the 
many particle case 0- 

The deterministic contribution, that reflects the behavior of the system for an infinite number of particles (N — > oo), 
can be described through a particle distribution p(x, t) that evolves according to the mean-field Fokker-Planck equation 
■ydtp(x,t) = [— a(t)d x F(x) + ksTd^] p(x,t). The net force per particle is a deterministic function f°°(t) = (F(x)) p = 

Jq dx p{x,t)F{x), with L the period of the ratchet potential. The net force is zero for the equilibrium distribution 
when the potential is on and also when it is off. We denote by f£° (t) with v — on, off the value of the deterministic 
part of the net force when the system has been evolving with the potential on or off respectively a time t after a 
switching. After a certain time t„ it can be approximately described by |j| 

f™{t)=C v e- x ^- r »\ (6) 

C v and Ay are constants that are obtained by fitting the net force obtained with the Fokker-Planck equation. In order 
to obtain (t) we evolve the equilibrium distribution for the off potential with the Fokker-Planck equation with the 
potential on, i.e., we assume that the system was close to the equilibrium state for the off potential before the off-on 
switching. We proceed analogously for f^(t). 

On the other hand, the amplitude of the fluctuations of the net force / can be estimated as Q 



This simple result is a good estimation of the amplitude of the fluctuations for potentials with characteristic height 
Vq and asymmetry a. For example, the potential 



[2ttx\ 1 /47raA 
Sin {—) + 2 S[n {—) 



(8) 



that we have used for the figures of this article, has characteristic height Vo and characteristic asymmetry a = 1/3 
(where ah is defined as the minimum distance between a minimum and a maximum of the potential, with L being 
the period of the potential). 

We have already provided estimations for both the deterministic part of the net force per particle and the amplitude 
of its fluctuations. This will allow us to calculate the average current. 

First, we compute the characteristic times during which the potential remains on, t OIL , and off, i ff. In the threshold 
control protocol the switching happens when the force crosses the threshold value with the appropriate slope [see 
Eq. ijlj]. When the threshold is crossed the equality u v = fu(t v ) is satisfied, with f u (t) the value of the net force a 
time t after a switching. Therefore, using the quasideterministic approximation JSJ wc obtain for the characteristic 
times 

|/~(**)|-£ = |u„|. (9) 
Using Eq. © we get the following explicit equations for the characteristic times 

\ v \U V \ + 2. 

with S given by Eq. Q. Moreover, Eq. (jSJ implies that this approximation is valid for t v > t„, where r„ are the 
transient times for each dynamics (afterwards, Eq. © is a good approximation). This implies Iw^l + X <C \C V \, 
that can be expressed as \u u \ + S <C max t \f£°(t)\ by using | | ~ max t \f£°(t)\. As S ~ 1/VJV, we see that this 
approximation is valid for small thresholds and large number of particles. 

We now compute the average displacement of the center-of-mass during an on-off period. Note that the center-of- 
mass moves only when the potential is 'on', because when it is 'off' the dynamics is purely diffusive. Therefore, as 
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the center-of-mass position is x cm = ^^Xi/N, its average displacement during an on-off cycle in the many particle 
case is given by using the evolution equations Q as 

Ax cm (t on ) = - [ "* f™(t) dt. (11) 
7 Jo 

The integration of the late time expression JSJl with v = on suggests a functional form 

Ax cm (ion) = Axon (l - e^"^") . (12) 

This functional form fits well the function Air cm (i on ) obtained from the numerical integration of the Fokker-Planck 
equation, and this fit is used to determine Ax on and Ai on - We have seen that the inclusion of the characteristic 
time At on improves the analytical results obtained in Ref. 7] (there it was assumed Ax cm (t on ) = Ax on ). This better 
estimation of the average displacement improves the results for the intermediate regime of not-so-large number of 
particles. Furthermore, the whole expression (|12|) is also necessary to improve the results for nonzero thresholds. 
When thresholds are enlarged the frequency of switching increases and therefore the times t on decrease. This implies 
a shorter displacement, as Eq. (fT2*)l predicts. 

The previous results allows us to give an approximate expression for the average center-of-mass velocity in the 
stationary regime, 

,. . .. x cm (t) - x cm (0) Ax oa / _ t /At \ Ax on [l-A(u on + EY/^ M ^] 

(arcm)st = hm = - II - e on/ J = j—— — — 5—— — — , (13) 

t^ 00 t t on +t oS \ / B — m(it on + S) — j^— ln(|ii ff | + S) 

with £ given by Eq. l(7)l. and A and B given by 

A = e -r,./Aw C 4/(A„Ag t B = Tm + ToS+ -L i n Con + J- In |C of f| • 

Aon A fF 

The final expression in Eq. I|13[l shows the explicit dependence on the thresholds u on , u g, and on the amplitude of 
the force fluctuations S; all the other parameters are determined by the dynamics for an infinite number of particles 
with zero thresholds. Eq. (|13fl has been obtained in the quasideterministic approximation and therefore is valid when 
the number of particles N is large and the thresholds are small as discussed after Eq. I|10(l . We have verified that 
it gives good estimations inside its regime of validity. In particular, for zero thresholds Eq. (|13fl is better than the 
formula obtained in Ref. Q thanks to the introduction of the characteristic time Ai on . (The formula in Ref. @ is 
recovered for u on = u g = and At on — 0.) 

Figs. Hl3*l compare the predictions of the quasideterministic approximation, Eq. I|13(l . with the numerical results for 
the threshold control protocol applied with the potential JSJ and Vq = 5ksT. For this potential the fit to the Fokker- 
Planck evolution gives C on = 0.67k B T/L, r on = 0.058L 2 /L>, A on = 28D/L 2 , C* off = -0.74fc s T/L, t oS = 0.037L 2 /D, 
A of f = 39D/L 2 , and Ax on = 0.08L, At on = 0.05L 2 /D. 

In Fig.^we plot the current as a function of the threshold u on (with u k = —u on ) comparing the quasideterministic 
approximation (|13|) and the numerical results obtained from the Langevin evolution equations We see that the 
quasideterministic approximation gives a good estimation of the current. However, it fails to predict the minimum 
located at low threshold values. This minimum is caused by a secondary effect that has not been accounted in 
the deduction of the analytic formula. This secondary effect is due to the fact that nonzero thresholds have the 
disadvantage of not being instantly optimal, because they imply switching on the potential when the force is still 
negative and switching off the potential when the force is still positive. In addition, for very small thresholds the 
switchings are not induced much earlier than they would be with zero thresholds due to the force fluctuations. Thus, 
there is a minimum located at thresholds of order l/v^/V, the magnitude of the force fluctuations. For larger threshold 
this secondary effect of the thresholds is overcompensated by the main effect of avoiding the undesired trapping of 
the dynamics. This main effect allows to have similar average displacements of the particles in a shorter on-off cycle 
time. Therefore, larger thresholds increase the average center-of-mass velocity. 

Figs. [21 and |3 compare analytic and numerical results for the current as a function of the number of particles for fixed 
nonzero thresholds: Fig.|21for u on = — u g = Q.lksT/L and Fig. for u on = 0.6ksT/L and u g — —OAksT/L (which 
are the optimal values for an infinite number of particles). In Fig. [3 we see that the quasideterministic approximation 
gives a good estimation for large number of particles. In Fig. the estimate is more rough due to the fact that the 
thresholds do not strictly verify the validity condition of the quasideterministic approximation (\u^\ + X <C |C„|). 
Another interesting result we have found is that for fixed nonzero thresholds the average velocity as a function of 
N tends to a constant asymptotic value for large number of particles, as Eq. (|13|) predicts. For an infinite number 
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FIG. 1: Average of the center-of-mass velocity (x cm )st as a function of the threshold u on for numbers 
of particles N = 10 5 ,10 6 and the limit N — ► oo for the potential © with V a — bk B T. An- 
alytical quasideterministic approximation (lines) and numerical results from Langevin equations i ffl l 
(points with error bars). We have taken u Q g = — u on . (Units: L = 1, D = 1, fcgT = 1.) 
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FIG. 2: Average of the center-of-mass velocity (i C m)st as a function of the number of particles N for the potential (jBJ 
with Vq — 5kgT and for thresholds u on = 0.1 and u s = —0.1. The simulations results obtained solving numerically 
the Langevin equations (points with error bars) are compared with the quasideterministic approximation for 
large N [Eq. l|13fl ] and the pure stochastic approximation for small N [Eq. i|19|) ]. The dotted horizontal straight 
line corresponds to the periodic switching protocol with optimal periods. (Units: L = 1, D = 1, kgT = 1.) 

of particles the force fluctuations vanishes, thus, this asymptotic value is given by Eq. (|13(l evaluated at S = 0. See 
Figs. 0and El 

The optimal threshold protocol gives the same current or better than the optimal periodic control @ (Fig. |3J). 
In particular, for an infinite number of particles the force fluctuations become negligible and the threshold control 
becomes equivalent to a periodic switching. The relation between the thresholds and the periods @ 

u v = f2°(T v ) (14) 

is obtained here as the limit N — > oo, i.e. £ = 0, of Eq. @. This relation permits to get the optimal thresholds for an 
infinite number of particles from the optimal periods just using the functions f^(t) and f^(t) obtained numerically 
from the Fokker-Planck equation. This avoids the need of integrating numerically N coupled Langevin equations for 
large values of TV. We have numerically checked that the expression (|14(l gives the optimal thresholds (see Sec. IIII Bl 
and Fig. 0|. 
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FIG. 3: Average of the center-of-mass velocity (i cm ) s t as a function of the number of particles N for 
the potential JSJ with Vb = 5fcsT and for thresholds u on = 0.6 and u s = —0.4 (optimal val- 
ues for N — > oo). The simulations results obtained solving numerically the Langevin equations 
(points with error bars) are compared with the quasideterministic approximation for large N [Eq. 1)13(1] 
and the pure stochastic approximation for small N [Eq. (|19|l ]. The dotted horizontal straight line corre- 
sponds to the periodic switching protocol with optimal periods. (Units: L = 1, D = 1, fc^T = 1.) 



2. Few particles: pure stochastic approximation 

When we have few particles the situation is the opposite to that considered in the previous section and the net 
force has nearly a pure stochastic behavior. A binomial distribution is found for the net force probability distribution, 
p(/), in Ref. |7j under the approximations that the position of the particles are statistically independent and that the 
probability of finding a particle in the interval [0, aL] is a. For simplicity this binomial distribution for the net force 
can be approximated by a Gaussian distribution 

*fl - 7m (15) 

with £ the amplitude of the fluctuations of the net force, that is given by Eq. J7J). Neglecting the time correlations in 
the net force, the average center-of-mass velocity for the threshold protocol [Eq. 10}] is given by 

(icm)st = - ^ Mf) df + ~ fP+(f) d f > (16) 

7 Ju aa 7 -'"off 

with p+{f) the probability of having a net force / and a non-negative value of / exp [p+(f) ~ p(/)/2]- This implies 
that, in the validity range of this small N approximation [S > max t |/°°(i)|], the current is a decreasing function of 
the threshold Won t CcLH be easily proven comparing the results for u' on and M on with < u' on < u on . Eq. gives 

1 [ u ° n 

(x cm )st{u' on ) - (x cm ) st (u on ) = - / /p_(/)c(f, (17) 

7 J u' on 

with p-(f) = p(f) — P+(f) > 0. Thus, the last term in the previous expression is non- negative implying 

(icro)st(^on) ~ (icm)st(Uon) > 0. (18) 

Analogously, it can be shown that for > u' oS > u oB we have (x cm ) st (u' oS ) - (i cm } st (u oS ) = (1/7) f™°® (- f)p+(f)df > 
0. This shows that the average center-of-mass velocity is a decreasing function for increasing modulus of the thresholds. 
Therefore, for small N we get the maximum current for zero thresholds. 

For small thresholds we have found an approximate explicit analytical expression for the current. If u g — — u on 
the contribution of the second integral in Eq. 116|) is generally small, because it is the integration of a nearly odd 



FIG. 4: Average of the center-of-mass velocity (x cm ) s t as a function of the threshold u on for 
N = 2,5, and 10 particles for the potential JSJ with Vq = 5ksT. Analytical pure stochastic ap- 
proximation l|19|l (lines) and numerical results from Langevin equations (points with error bars) 
are compared. We have taken u Q g = —u on . (Units: L = 1, D = 1, kgT = 1.) 

function in a nearly symmetric interval around zero. On the other hand, the contribution of the first integral is greater 
provided the thresholds are small (u on < £). Then, neglecting the second integral we obtain 

S "on 

(icm)st - — 7S= e ~^- ( 19 ) 
7V 2n 

(Note that for u on = we recover the zero threshold result found in Ref. 0.) This expression, Eq. I|19|) . gives good 
predictions when we have few particles and small thresholds. In particular, we show in Figs. 12141 that it correctly 
predicts the threshold and particle number dependence of the current, even for u on ~ £ ~ 3.4 when N = 10 (Fig. 0}. 

B. General thresholds 

In the previous subsection we have studied the threshold protocol when the moduli of the thresholds are small, 
obtaining approximate analytical expressions for the current. In contrast, in this subsection we study the threshold 
protocol for general thresholds (that are in general beyond the applicability range of the previous analytical expres- 
sions) . This study is done performing numerical simulations of the Langevin equation of the threshold protocol for 
general values of the thresholds. 

1. U ff = —u on 

Let us discuss first the results for thresholds that are related by u ff — — "on- 

In the few particle case, when the thresholds are small the current decreases exponentially with the square of the 
threshold as we have already seen [see Eq. (|19|) ]. However, as the rate of the exponential is small, we nearly have a 
plateau near the maximum at zero thresholds, as shown in Figs.0]and|SJ On the other hand, for very large thresholds 
Eq. (|19fl is no longer valid and the current decreases faster than the exponential. Note that the current continues to 
be a decreasing function, as predicted by Eq. I|18fl (valid for any threshold values in the few particle case). See Figs. 01 
and [3] 

In contrast, in the many particle case the maximum of the current is no longer at zero thresholds, but at a finite 
value. As we have explained before, the introduction of thresholds has the advantage of inducing earlier switchings. 
This avoids the undesired trapping that otherwise is present for large N implying low current values. The presence of 
thresholds allows to have similar average displacements of the particles in a shorter on-off cycle time, and therefore 
increases the average center-of-mass velocity. However, if the thresholds are too large the losses in the displacement 
become more important than the gains of shortening the on-off cycle time. Therefore, the current has a maximum 
located at a finite value of the thresholds in the many particle case (Fig. |SJl . (The tiny minimum in the small threshold 
region is related to another effect: the disadvantages of choosing a not instantly optimal protocol. For a more detailed 
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FIG. 5: Average of the center-of-mass velocity (x C m)st as a function of the threshold u on with 
UqB — —u n for various N. The lines correspond to the numerical solution of the Langevin equa- 
tions (HJ for the potential © with V = 5k B T. (Units: L = 1, D = 1, k B T = 1.) 

explanation see Subsec. Till A II ) 

Another important result in the many particle case is that the maximum obtained for the current as a function of the 
threshold magnitude is quite flat and nearly independent of the number of particles. See Fig- EI 

In summary, in the many particle case the current has a maximum for nonzero thresholds whose position is nearly 
independent of the number of particles. On the other hand, for few particles the current is maximum for zero 
thresholds. However, in the few particle case the current is nearly the same up to thresholds of the order of the 
thresholds that give the maximum for the many particle case (see Figs. 0] and |SJ. This has an important implication: 
the optimal thresholds values for the many particle case give currents close to the maximum for any number of 
particles. 

2. u ff / -u on 

The study of the current for completely general thresholds u on and u s (without restrictions) reveals that the 
behavior is analogous to that described previously. In fact, the optimal thresholds for large number of particles are 
located not far from the line u on — — u ff, and these thresholds give currents close to the maximum for any number 
of particles. (See Figs. and 0) 

As we have already commented in the previous section, for an infinite number of particles the force fluctuations 
becomes negligible and the threshold protocol becomes equivalent to a periodic switching. This implies the relation 114|) 
between the optimal periods T on and T Q g, and the optimal thresholds M on and ii ff, that we have numerically checked 
(see Fig. UJ. Therefore, these relations permit to obtain the optimal thresholds for an infinite number of particles 
from the optimal periods, just using the functions f^(t) an d /<SfW obtained numerically from the Fokker-Planck 
equation. These thresholds give good results for large number of particles. Moreover, it is important to note that these 
thresholds values also give currents close to the maximum in the few particles case due to the smooth dependence for 
small thresholds (see Figs. I5I7|) . 

In particular, we have seen that for the potential JSJ) with Vq = bk B T the optimal switching periods are ap- 
proximately T on = O.Q6L 2 /D and T Q g — 0.05L 2 /D. Therefore, with just a Fokker-Planck simulation for the poten- 
tial we have found that a good estimation of the optimal thresholds is given by w on = /^(Ton) — 0.6k B T/L and 
Uoft — ftff(%ff) = —0Ak B T/L, in good agreement with Fig.0 

IV. CONCLUSIONS 

In this paper we have analyzed the threshold control protocol for a collective flashing ratchet. We have studied the 
threshold dependence of the current in this closed-loop control protocol. The quasideterministic (for many particles) 
approximation 7] has been improved through the introduction of an additional characteristic time giving better results 
for not-so-many particles. Both the quasideterministic and the stochastic (for few particles) approximations Q have 
been applied to the threshold control protocol. This has led to analytical expressions for large and small number of 
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FIG. 6: Thresholds dependence of the average of the center-of-mass velocity (i cm )st for N = 10 4 par- 
ticles in the potential JSJ with Vo = hksT. The grid has been obtained integrating numerically 
Langevin equations Q for different thresholds u on and w Q ff (Units: L = 1, D = 1, fc^T = 1.) 




FIG. 7: Thresholds contour lines corresponding to the value of the average of the center-of-mass velocity (i C m)st 
5% below its maximum for N = 10 2 , N = 10 3 and N = 10 4 particles in the potential © with V = 5k B T. 
The contour line for N = 10 5 is already very similar to that for N = 10 4 . The point corresponds to the optimal 
thresholds for N — ► oo obtained from the optimal periods using Eq. I|14|). (Units: L = 1, Z) = 1, /j^T = 1.) 



particles. We have computed numerically the current dependence on the thresholds and on the number of particles 
obtaining a good agreement between analytical and numerical results in the validity range of our assumptions. We 
have also compared these results with the optimal periodic switching protocol. 

We have seen that for many particles the current has a maximum for nonzero thresholds whose position is nearly 
independent of the number of particles. On the other hand, for few particles we have demonstrated that the current 
increases as thresholds moduli decrease, so the maximum current is reached at zero thresholds. However, the current 
is nearly the same up to thresholds of the order of the optimal thresholds for the many particle case. This implies 
that the optimal thresholds values for the many particle case give currents close to the maximum for any number 
of particles. The optimal thresholds for an infinite number of particles can be obtained from the optimal periods 
of the periodic protocol just solving the Fokker-Planck equation in two particular cases (potential 'on' and 'off', see 
Section fill A() . Therefore, we can get a good estimation of the optimal thresholds for many particles, that also gives 
currents close to the optimal for any number of particles as we have shown. 

The closed-loop threshold control gives the same current as the optimal protocols for the one particle case and for 
an infinite number of particles, and it gives high currents in between. However, obtaining the best protocol for the 
maximization of the current is still an open question. 

In this work, and in previous ones 0, @j we have seen that thanks to the information about the fluctuations 
obtained through the feedback, the performance of the system can be increased. This increase of the performance 
has thermodynamical limitations that have been studied in a general context for the efficiency 10]. We are now 
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working in order to get a deeper understanding of this interplay between the information and the increase of the 
performance |ll|. 
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